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Hall viscosity from effective field theory 
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For two-dimensional non-dissipative fluids with broken parity, we show via effective field the- 
ory methods that the infrared dynamics generically exhibit Hall viscosity — a conservative form of 
viscosity compatible with two-dimensional isotropy. The equality between the Hall viscosity coef- 
ficient and the ground state's intrinsic angular momentum density follows straightforwardly from 
their descending from the same Lagrangian term of the low-energy effective action. We show that 
for such fluids sound waves are not purely longitudinal, but acquire an elliptical polarization, with 
transverse-to-longitudinal aspect ratio proportional to frequency. Our analysis is fully relativistic, 
thus providing a natural description of (2+1) dimensional relativistic fluids with broken parity. 



Introduction. At large distances and long times, the be- 
havior of any fluid can be described by the equations of 
fluid dynamics — the continuity equation and the Navier- 
Stokes equation. These equations are in essense the equa- 
tions of mass and momentum conservation, 
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coupled with an equation expressing the stress tensor T %3 
via the p and v l . The latter can be written down based 
on symmetry considerations and involves shear and bulk 
viscosities. 

In two dimensional fluids, it has been noticed some 
time ago [l|, Q that new terms can be added into the 
hydrodynamic equation if one relaxes the condition of 
parity and time reversal invariance. Namely, one can 
include into the stress tensor the Hall viscosity, which is 
odd under these discrete symmetry, and enters the hydro- 
dynamic equations at the same order as the conventional 
viscosities. By nature, the Hall viscosity is dissipation- 
less. 

The Hall viscosity has been investigated mostly in 
gapped systems, where it is computed by subjecting the 
system under consideration to a shear metric perturba- 
tion which is slowly changing with time. The Hall vis- 
cosity has been shown to be related to a Berry phase. It 
was found that the Hall viscosity is proportional to the 
density of intrinsic angular momentum: t}a = \sn where 
n is the particle number density and s is the average spin 
per particle (so sn is the area density of intrinsic angular 
momentum). 

In this paper we construct a theory of a two- 
dimensional compressible fluid with broken spatial par- 
ity. We construct such a theory by extending the action 



of a perfect fluid to include parity-breaking effects. As 
the theory is based on an action principle, it does not 
contain dissipation effects. We find that, when the lead- 
ing parity-breaking term is included into the Lagrangian, 
the equations of motion describe a fluid with a nonzero 
intrinsic angular momentum density. In this fluid, in 
order to have a symmetric stress tensor, in general the 
momentum density must be defined to contain a term 
proportional to the derivatives of the intrinsic angular 
momentum density, in addition to the pv term. 

However, we find that in one particular case, namely 
when the intrinsic angular momentum density is propor- 
tional to the particle number density (i.e., the intrin- 
sic angular momentum per particle is constant), there is 
an alternative formulation of the hydrodynamic theory 
where the momentum density remains pv, but the stress 
tensor now contains the Hall viscosity. The value of the 
Hall viscosity is exactly half of the intrinsic angular mo- 
mentum density. 

Thus, in our compressible fluid, the Hall viscosity can 
be thought of as an alternative description of a two- 
dimensional fluid with a constant intrinsic angular mo- 
mentum per particle. 

The setup. We will use the field-theoretical description of 
fluids introduced in Ref. [3|, which is based on previous 
work on the field theory of solids 4, 5] and supersolids [61 . 
We refer the reader to Ref. [1] for details (see also Refs. [3, 
8]). The starting point is the parameterization of a fluid's 
configuration space by giving at time t the comoving (or 
"Lagrangian" ) coordinates (j) 1 as functions of the volume 
elements' physical (or "Eulerian") positions x: 



= <p T (x,t) 



1 = 1, 
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where d is the number of spatial dimensions. This is com- 
pletely equivalent to the inverse point of view where one 
gives the physical positions x as functions of the comov- 
ing ones, x = x(4> r ,t) (this parametrization is adopted, 
e.g., in Ref. Q]). The advantage of using §5§ is that cj) 1 can 
be treated as scalar fields in a field theory with spacetime 
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symmetries. We can thus construct local Lagrangians 
with the desired symmetries following the usual rules of 
(effective) field theory. 

In a solid, the natural choice for the comoving coordi- 
nates 4* 1 is to point along the crystal axes. For a fluid, 
there is an arbitrariness in assigning comoving coordi- 
nates to the individual volume elements. If the fluid is 
at rest in the infinite past at constant density po, then 
the most convenient choice is to identify the comoving 
coordinates with the physical ones at t —> — oo 

(j) 1 = x 1 . (4) 

With this choice of coordinates, the internal symmetries 
the dynamics must obey are Q 



b 1 -> 6 1 



a 

R T j<i> J 



1=1,. ..,d 

R E SO(d) 

A* 9 *' 1 

det 7u7 = 1 



(5) 
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Equations. ([3]) and ([5]) correspond to the physical homo- 
geneity and isotropy of the fluid's internal space. Eq. ([7} 
is what distinguishes a fluid from an isotropic solid: dis- 
placing volume elements without compressing or dilating 
the fluid anywhere does not cost any energy. 

Our main goal will be to construct parity-breaking hy- 
drodynamics for non-relativistic fluids, but we find it con- 
venient to keep our analysis fully relativistic and to take 
the non-relativistic limit only when needed. Our symme- 
tries are therefore eqs. ([SHU) plus the (d+ l)-dimensional 
Poincare group, under which the 0*'s behave as scalars. 
We will use the (—,+,...,+) signature for the metric. A 
fundamental object that we will use extensively is 
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(We define the [d + l)-dimensional e tensor by e 01 " = 
+1.) It is a vector under Poincare, and is invariant under 
our internal symmetries ([5]-[7]). We will denote its norm 
by b, 



b= tJ-J^J* 1 = VdetB 17 



B 



u 



d^d^ J , (9) 



which gives a measure of the compression level of the 
fluid. The velocity field u M of a fluid configuration is 
defined as a unit timelike vector aligned with 



= bit* 1 . 

The current J M is identically conserved: 
d^J^ = (identically). 



(10) 



(11) 



Equations (|10|). (fTTj) invite the interpretation of b as the 
number density of fluid elements, or of fluid points, whose 
conservation law should indeed be an identity. 



The Lagrangian and the stress-energy tensor. According 
to standard effective field theory logic, at low energies 
and momenta we should organize the dynamics as an ex- 
pansion in derivatives. The low-energy Lagrangian thus 
takes the form Q 



C = F(b) + higher derivatives , 



(12) 



where F is a generic function, to be determined from the 
equation of state. More precisely, from the lowest deriva- 
tive part of the Lagrangian, F(b), one gets the standard 
perfect fluid stress energy tensor, 

Tf = (p + p)u a ^+ PV ^, (13) 

with energy density and pressure given by Q 

p = -F(b) , p = F{b) - F'{b)b . (14) 

For a barotropic fluid the equation of state is an algebraic 
relation between p and p: p — p(p). This translates, via 
Eq. (fT4"|) , into a differential equation for F. 

The higher derivative terms generically involve at least 
two more derivatives, on top of those already present in 
objects like J M and b, which involve one derivative per 
field. Indeed, the shift symmetry (J5J) forces each field (j) 1 
to be acted upon by at least one derivative. So, at lowest 
order in derivatives the only object that is invariant under 
the internal symmetries and that transforms covariantly 
under Poincare is J M itself, along with functions thereof 
like b. If we add only one more derivative, the general 
structure we expect is a generic function of 6, times a 
derivative on J M , times extra powers of J M , 



AC = f(b) dJJ n 



(15) 



with suitable contractions of the Lorentz indices. No- 
tice that a derivative acting on a function of J M is still 
rewritable in the above form, via the chain rule. The 
derivative cannot be contracted with the J it acts on, 
because of ([TT]) . Moreover, all contractions between two 
J's both belonging to the J" piece just redefine f(b), 
since J^J^ = —b 2 . Finally, the J n piece cannot have 
more than one J contracted with an e tensor, because 
of symmetry. So, for generic dimensionality d, at the 
one-derivative level the only possibility is 



ac = f(b) d^j v = d^b 2 j» 

= d^g(b) J» , 



(16) 



13j . One is 



which vanishes upon integrating by parts 
thus led to consider terms with more derivatives. 

However in 2 + 1 dimensions we have one more possi- 
bility: 



AC = f(b) e^ p J^d„J p , 



(17) 



This is invariant under all our symmetries, and should 
generically be there in two-dimensional fluids where par- 
ity is broken. Its contribution to the stress tensor will 
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involve a single derivative acting on physical quantities, 
like the velocity field for instance. In this sense it is of 
the same order as viscosity — it appears at the same or- 
der in the derivative expansion. However, coming from 
a Lagrangian term, this will be a conservative (i.e., non- 
dissipative) form of viscosity. Our claim, which we are 
now going to prove, is that (I17|) is the Lagrangian de- 
scription of Hall viscosity. 

To compute the stress-energy tensor we vary AC with 
respect to the metric. After a straightforward computa- 
tion we get 

AT* ee 2 ^ = AC {rf* + »V) 

og a p a log b 

+ Afe u ^ a jVd v J ll + 2d v fM a jVj ll (18) 

where we made use of the identity 

BJj d a cj> I d fS (j> J = n afs + u a u p . (19) 

At the one-derivative level, the full stress energy tensor 
is 

T a/3 = jaP + AT a/3 ^ (20) 

where Tq^ is the perfect fluid part, eq. (fT3|) . 

Intrinsic angular momentum and Hall viscosity. A cru- 
cial fact is that the otherwise arbitrary coefficient fib) 
turns out to be related to the angular momentum density 
for static configurations. The total angular momentum 
is defined as 

L = e tJ Jd 2 xx l T 0j (21) 

(it is a two-dimensional scalar.) For vanishing fluid ve- 
locity, T° l reduces to 

T M = AT oi = -e y ' dj (fb 2 ) . (22) 

Upon plugging this into the expression for L, integrat- 
ing by parts, and discarding the boundary term which 
vanishes for a finite-size fluid (if we take the boundary 
outside the fluid), we get an angular momentum surface 
density 

l^ d ± = -2fb 2 . (23) 

The fact that the angular momentum density does not 
vanish whhen the fluid is at rest is not surprising — parity 
breaking allows it. This situation occurs, for example, in 
the A-phase of superfluid helium-3 [ll[ . 

We now take the nonrelativistic limit. In this limit, 
b is the mass density /?, so we will use b and p inter- 
changeably [H} • At lowest order in the fluid velocity, our 
corrections read 

AT 00 = O(dlv) (24) 

AT 0t = dj£ + 0{d£v 2 ) (25) 

AT 13 = -(£-±b£') (e kl d k vi) 5 ij 

+ i( £ , VW + mj)+0(9fo 3 ), (26) 



where we parameterized everything in terms of the an- 
gular momentum density £(b) = —f(b)b 2 , £' stands for 
the derivative of £ with respect to b, and the d inside the 
C(...)'s denotes schematically a spatial gradient. One 
can check that all 0(. . .) are suppressed in the nonrela- 
tivistic limit and can be dropped. 

First, we see that AT y does not resemble the stress 
tensor associated with Hall viscosity. Its traceless part 
contains terms proportional to the gradient of the angular 
momentum density £, but not gradients of the velocity. 
Moreover, the non-relativistic momentum density is not 
the naive pv l , being corrected by a term proportional to 
the gradient of the angular momentum, 

T 0t = pv i + \e 13 dj£ . (27) 

That the momentum density involves the gradient of 
the angular momentum density is nothing new [III [l2j . 
Galilean invariance then dictates that the stress tensor 
must contains terms of the form vd£, as in Eq. (I2"r?l) (in 
contrast, the term proportional to e kl dkVi in Eq. ([2U)) is 
not dictated by Galilean invariance). 

The extra term in the momentum density does not 
contribute to the total momentum J d 2 x T° l , and has zero 
divergence: diAT° l — 0. For this reason, the continuity 
equation remains (fTJ), and the momentum conservation 
can still be written as Eq. ([2]), but with a new stress 
tensor, 

fit (pv l ) + d, (T IJ + E lJ ) = , E tf = - le ij £' d k {b v k ) 

(28) 

£ y has been defined so that <9;AT° J = fi,S y . The mod- 
ified stress tensor is not symmetric. It should be ex- 
pected since the symmetry of this stress tensor would 
lead to the conservation of the naive angular momentum 
j dx p e 8J , but we know that the conserved angular 
momentum is actually the sum of the naive angular mo- 
mentum and the "spin," 

Jdx{pe ij x l v J +£) (29) 

In one particular case, however, we should be able to 
symmetrize T, lJ via the Belinfante trick: when the an- 
gular momentum density is proportional to the particle 
number density: 

£ = Xb , A = const . (30) 

In this case the total "spin" is proportional to the to- 
tal particle number and is conserved by itself. This is a 
non-trivial assumption, which is obeyed for instance by 
fluids where the bulk of the angular momentum is car- 
ried by individual spins or by bound states, so that in 
first approximation it scales linearly with the number of 
particles in the system. Under this assumption we can 
symmetrize U* 3 via the addition of the total divergence of 
an antisymmetric tensor, c^A 1 ^, which does not affect 
momentum conservation. We have to choose 

A* fc l = iAfep' v k - e lk v j - e jk i/] , (31) 
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from which we get 

S lJ + fe A* fc ] = -i\[ e ik d k (bv j ) (32) 

which is symmetric, as desired. Putting everything to- 
gether we get 



f 00 = p, 
f m = pv i , 

f ij = P 5 ij +pv i v j + AT ij 

with 

AT ij = AT ij + S lJ + d k A i[]k] 

= -|A& [(e kl d kVl ) 5* + {e lk d k v' + i j)] 



(33) 
(34) 
(35) 



(36) 



Notice that all derivatives of the number density b can- 
celed out, leaving us with derivatives of the velocity field 
only. The form of AT y matches precisely Hall viscos- 
ity [H, • It can be rewritten as 

ATij = -m(^ikSji + ej k S t i)V k i , Vm = \{d k vi + d t v k ) 

(37) 

with 



m = §a& = ¥ ■ 



(38) 



The Hall viscosity coefficient is thus half the angular 
momentum density. This relationship between the Hall 
viscosity and the angular momentum density has been 
shown to occur in gapped system using adiabatic argu- 
ments [l^ . In our description of a compressible fluid, it 
emerges as a straightforward consequence of the simplic- 
ity of the low-energy effective action. At next-to-lowest 
order in the derivative expansion, all observables derive 
from a single Lagrangian term, eq. (|17p . Moreover, our 
Lagrangian provides a relativistic generalization of Hall 
viscosity. 

Sound wave propagation. We now investigate sound-wave 
propagation in a compressible fluid with broken parity. 
Sound waves are obtained from our field theoretical de- 
scription by expanding the Lagrangian at quadratic order 
in small perturbations. Explicitly, if we consider small 
deviations from the static, homogeneous configuration 
©, 

4> a = x a + ir a , (39) 

where 7r a is (up to a sign) the displacement, the leading- 
order Lagrangian, to second order in the phonon field tt, 
becomes 0, Q 



F(b) -> w [±t? 2 



(40) 



where wq = —F'(l) = (p + p)b=i is the background 
enthalpy density (which becomes the mass density in 
the nonrelativistic limit), and <? s = F"{\)/F'{\) — 
(dp/dp)b=i is the sound speed. As expected, only the 
longtidinal part of ir propagates. The transverse modes 



lack gradient energy — this is a direct consequence of our 
symmetry ([7]) — and do not propagate. 

We now restrict to d — 2 case, add our one-derivative 
correction, Eq. (|17j) . and expand it at second order in 7?. 
By using eq. (fT0|) and 

U = I +C(7T 2 ), U = -7?+(7T- V)7? + 0(7T 3 ), (41) 

b= l + V-7f + 0(7r 2 ), (42) 



we get 



AC 



e'd k ir k e ij diic j ] 



(43) 



where I and £' are evaluated at b = 1 , we used the univer- 
sal relation between / and £, eq. (|2"5)l . but we did not as- 
sume Eq. (|30[) . nor did we take the non- relativistic limit. 
This correction to the quadratic Lagrangian (|4"0")) induces 
a mixing between longitudinal and transverse modes, as 
can be immediately seen by decomposing the phonon field 
as 7? = 7?£ + ttt- To diagonalize the Lagrangian, one can 
introduce a new phonon field 7r', with longitudinal and 
transverse components if' L and tt' t , respectively, such that 



2wq cl 



1 
it'll 



(44) 
(45) 



One gets simply 



C 



wo\ 



\c 2 s (V -tt') 2 ] +higher derivatives , (46) 



where now the higher derivative terms involve at least two 
more derivatives, i.e. they are on an equal footing with 
higher-order corrections that we have been neglecting all 
along. It is thus consistent to discard them — and in fact 
it would be inconsistent not to. We see that at the order 
we are working the spectrum of perturbations is unal- 
tered: n' T does not propagate while if' L features a linear 
dispersion relation with propagation speed c s . However 
what is interesting is that this propagating mode does not 
look purely longitudinal once expressed in terms of the 
original field ir: it is elliptically polarized, alternating be- 
tween the longitudinal and the transverse direction, with 
an aspect ratio proportional to frequency, 



w 



2c1 



-I 



(47) 



We do not expect generically a cancellation between the 
two terms inside the bracket: one depends on the sound 
speed while the other does not, and, in particular, for 
a non-relativistic fluid the latter will be negligible with 
respect to the former. 

Conclusions. In this paper we have considered a field 
theory description of a (2+1) dimensional fluid with bro- 
ken parity. We have shown that the fluid naturally has 
a finite density of angular momentum, and that in the 
nonrelativistic limit it allows a description involving Hall 
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viscosity if the intrinsic angular momentum per particle 
is constant. The Hall viscosity is shown to be exactly half 
of the area density of the intrinsic angular momentum. 
We also investigated the propagation of sound wave in 
such a fluid. While we did not find modifications to the 
sound speed, we found that the sound wave is a super- 
position of transverse and longitudinal waves. 

Although in this paper we have taken the point of view 
that the parity-breaking effect is the next-to-leading or- 
der effect in the derivative expansion, one can regard the 
Lagrangian C + AC, without higher order terms, as a 
complete description of a hypothetical fluid, even when 



the AC is of the same order as C. For example, one can 
consider a fluid of particles with very high spin. The 
hydrodynamic equations that we derived do not suffer 
inconsistencies when truncated to the order that we have 
considered. 
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